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ABSTRACT 

A scheme of quantum electrodynamic (QED) background-free radiative emission of neutrino pair (RENP) 
is proposed in order to achieve precision determination of neutrino properties so far not accessible. The 
important point for the background rejection is the fact that the dispersion relation between wave vector 
along propagating direction in wave guide (and in a photonic-crystal type fiber) and frequency is modified 
by a discretized non-vanishing effective mass. This effective mass acts as a cutoff of allowed frequencies, 
and one may select the RENP photon energy region free of all macro-coherently amplified QED processes 
by choosing the cutoff larger than the mass of neutrinos. 
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1 Introduction 


With the advent of successful macro-coherent amplification (more than 10^®) of QED rare processes [2], 
namely the macro-coherent paired super-radiance (PSR) [3], the atomic project of neutrino mass spec¬ 
troscopy [1] has gained a new stage of potentiality to explore important neutrino properties yet to be 
measured and to ultimately detect the relic neutrino of temperature 1.9 K [5]. 

In the present work we address the problem of quantum electrodynamic (QED) backgrounds and propose 
a scheme of QED background-free RENP (Radiative Emission of Neutrino Pair). Usual higher order QED 
processes, when they occur spontaneously, are not at all serious backgrounds to macro-coherently amplified 
atomic neutrino pair emission (RENP) if experiments of the neutrino mass spectroscopy are designed with a 
repetition cycle sufficiently faster than the decay lifetime. (Even a repetition scheme slower than the decay 
rate is conceivable if the dead time is not too large.) The macro-coherently amplified QED process may 
however become a serious source of backgrounds, since their rates are much larger, as is made evident below. 
We shall term macro-coherently amplified QED process of order n as McQn for brevity. The case of n = 2 
corresponds to PSR. Since RENP process occurs with parity change, the main backgrounds are odd McQn. 

Our proposal for the QED background rejection is to use either wave guides [6] or some type of pho¬ 
tonic crystals [7] to host a target. A promising host is Bragg fiber consisting a hollow surrounded by two 
periodically arranged dielectrics of a cylindrical shape [S], M- Eor brevity we call these hosts as host guides 
in the present work. After we show below how QED backgrounds are rejected, we shall calculate spec¬ 
tral rates of RENP (stimulated single photon emission) |e) —)• l^f) +^0 + 1 ^ 1 ? in which no background of 
McQ3 |e) —)• \g) -|- 70 + 7172 exists. Rejection of McQ5 and so on is then automatically guaranteed. The 
background-free RENP rate is, to a good approximation, found to be a shifted (to the higher energy side) 
spectrum in free space. 

Parities of two states, |e), \g), for RENP are different. A good example of candidate de-excitation is from 
the Xe excited state of = 1“ of configuration 5j?®6s(8.437 eV) (the decay rate being ~ 300 MHz) to the 
ground state of O"'' of 5p®. 

We show a part of Feynman diagrams (despite of the use of the non-relativistic perturbation theory 
based on bound state electrons) in Fig([T]). RENP diagrams are based on the nuclear monopole contribution 
of |10] . Relevant Xe energy levels are shown in Fig([2]). 
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Figure 1: A part of Feynman diagrams for Xe ^Pi(8.437eV) RENP in the left and for McQ3 in the right. 
Dashed red line in the left is for Coulomb excitation between nucleus A and a valence electron. 
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Figure 2: Relevant Xe energy levels: |e) = 6 s^Pi , \g) = 5 p^Sq , |g) = 7s^Pi ( 6 s^P 2 is for the spin current 
contribution in Appendix). Solid blue lines are the route for RENP, while dashed red lines are for McQ3. 


The rest of this paper is organized as follows. In Section 2 the main background source of QED processes, 
McQ3, is discussed in detail and its rate is given along with some features of background events. Results 
are presented for the Xe de-excitation and the necessity of all QED background rejection is emphasized. In 
Section 3 we point out the possibility of the QED background rejection in targets hosted in wave guides 
or photonic crystals, and derive a condition for the background rejection in RENP. The idea is based on 
the difference of the very nature of freely propagating neutrino and restricted propagating light field in the 
hosted target. The background-free RENP spectrum is approximately found to be the shifted spectrum in 
free space, the shifted threshold being determined by the frequency cutoff in host guide. A few examples of 
the shifted spectrum due to the nuclear monopole contribution m are shown taking a few choices of the 
size of host guides. In the first Appendix some rudimentary facts on propagating modes in wave guides and 
photonic-crystal fibers are explained. The second Appendix presents the modified RENP spectrum from 
6 s^P 2 (8.315 eV) due to the spin current contribution [TT] . 

Throughout this paper, we shall use the natural unit of ^ = c = 1. Conversion factors are thus 
(leV)-^ = 1.97 X lO-^cm = x lO^'^Hz)-!. 

2 McQ3 rate 

We shall first demonstrate that without an experimental method of background rejection McQ3 process 
gives an enormous background rate when the process is macro-coherently amplified. Thus, it becomes a 
necessity to invent the rejection method when one wants to detect RENP. 
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Kinematics 


For a fixed trigger frequency wq (of less than tegl‘1-, with the level spacing, as a necessary condition 
for the macro-coherence), there is a unique relation between the detected photon energy and its emitted 
angle 0*: 


cos 0, 


= — 1 — 
Wo 


li^eg 2wo) 


i = 1,2. 


( 1 ) 


This relation may be derived from the 4-vector relation of a — /cq = —{ki + ^ 2 ) where a = (eeg,0) is the 
energy-momentum 4-vector relevant to atomic transition (the difference between the initial and the final 
states), and ki,i = 0,1,2 are those of three photons. Note that with the macro-coherence amplification 
both of the energy and the momentum are conserved in the atomic process [3]. The angular constraint 
I cos 9i\ <1 gives ^ — wq < Wj < ^ . 


Basic data: 6s^Pi(8.437 eV) de-excitation 

The initial metastable state |e) is an electron-hole system consisting of a valence electron of 6 s and a 
hole of 5p or a filled state 5p^, their quantum numbers being much like those of two valence electron system. 
Relevant energy levels used are 

|e) = 6s3Pi(8.437eV), \p) = 6p^P2{9.821eV) , ( 2 ) 

k) = |e), I 5 ) = 5p^5o(OeV). (3) 

We have used the notation of LS coupling scheme, although other coupling scheme may be more appropriate 
for Xe. Two decay rates or A-coefficients for McQ3 rate computation are known: 

A( 6 s — 5p) (of Ig) —)• I 5 )) = 2.81 X 10®Hz , A{6p —)• 6 s) (of \p) —)• \q)) = 2.24 x lO’^Hz . (4) 


McQ3 rate in free space 

We need two more atomic levels in addition to \e),\g), to induce McQ3: they are denoted by |p),|(i'). 
These levels are connected to each other by electric dipole transitions of moments, dep,dpg,dqg. Quantum 
mechanical amplitude of McQ3 for a single route via virtual levels \e) ^ \p) ^ \q) ^ \g) is in proportion to 

dep ■ Padpq ' Ejydqg • Eq , (5) 

where a,b,c are taken respectively from (0,1,2). The magnitudes of dipole moments may be replaced by 
A-coefficient divided by the third power of their level spacings, hence in principle measured data can be 
used. 

In the example of Xe de-excitation, three photon emission occurs in the third order of QED: 

(6s5p“)j=i ( 6 p 5 p“)j =2 (6s5p“)j=i -F 7 -k 7 (5p^)j=o + 7 + 7 + 7, ( 6 ) 

where 5p~ denotes a hole in bp state. The final state denoted by bp^ is actually the completely filled six 
bp states 5p®. Three emitted photons here are taken from 70 , 71 , 72 ; hence there are 3! = 6 diagramatic 
contributions. 
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The differential McQ3 spectrum for a single detected photon of energy uj = oji summed over polarizations 
is calculated as [12] 
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( 8 ) 


^pg ^0 ^ ^ ^qg ^0 

The spectrum shapes are illustrated for two choices of wo taken close to eeg/2; at eeg/2 and its 5% reduced 
energy in Fig([3|). Despite of a non-trivial matrix element oc T 3 (w), the spectrum shape is symmetric around 
the point w = (cgg — wo)/2. This McQ3 rate is expected much larger than the RENT rate given below. 


Xe McQ3 spectrum 



Figure 3: Xe McQ3 spectrum rate in arbitrary unit. The actual rate values integrated over photon energies 
are (5,6) x lO^^Hz (n/10^‘^cm“^)^(R/cm^)(?73(t)/10“^) for these two cases of the trigger energy wo: wq = 
eeg/2 in solid black and 5% reduced value in dashed red. 


It is important to mention that McQn processes, in particular, the McQ3 case has an interesting poten¬ 
tiality of creating a quantum entangled pair of photons under the influence of the trigger, although this is 
not the appropriate place to discuss the subject in detail 

McQ3 in wave guide 

Although photonic-crystal fibers such as the Bragg fiber [8] are more promising due to their smaller 
transverse sizes and other reasons, we shall explain our ideas using the wave guide because concepts are 
simpler there. A good comparison of these two materials is given in [7|, in particular in its chapter 9. 

As explained in Appendix, propagation of light field in wave guide is restricted: compared to propagation 
in free space, allowed modes in wave guides have discretized transverse wave vectors, k±, and their smallest 
cutoff. 
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Due to this restriction, the continuous momentum integration over 3-vector k, is replaced by a semi¬ 
discrete sum. In the square wave guide of size d these discrete modes are given by wave vectors in two 
dimensions, 

k_L = ^{ni,n2), rii = 0,1,2,--- , -F n| > 1. (9) 

The discrete sum over transverse momentum components may however be approximated, , when deeg S> 
1, by a formula of the continuous momentum integration, however with the important difference in the 
presence of the transverse cutoff at kc = Tr/d. 


3 Background-free RENP rate 


Effective cutoff in host guides 

The dispersion relation between the energy oj and the longitudinal momentum A:|| along the propagation 
axis is given by w = in host guides. In the case of square wave guide of size d the effective 

mass M is discretized as M = \Jni + ri^ kc ,kc = n/d with rii integers, as explained in Appendix. A similar 
effective mass exists in Bragg fibers along with the band gap [8]. It is important to note that the effective 
mass M cannot vanish. 

A derivation of the condition for the QED background rejection in RENP is given for the general case 
of McQn (n > 3). We use the Lorentz invariance of Aq — in the temporal and the spatial components of 
4-vector A, (Aq, A^), where z-axis is taken along the propagation direction. We use the energy-momentum 
conservation: a — ko = —k — to give a scalar invariant relation among 2-vectors, 

n—1 

{a - kof = {k + '^ kif . (10) 

i=2 

The invariant left hand side (LHS) quantity may be evaluated in the laboratory frame, to give e^g — 2ujQeeg + 
Mq , while the right hand side (RHS) quantity is calculated in the CM system of n — 1 longitudinal momenta, 
to give (w -h Hence, 


n—1 

^Ig - 2a;oeeg -|- Mq = (w -|- ^ \/ ^‘i,z + ■ (11) 

i=2 

For the largest cuq value of allowed region, one takes the smallest of RHS, namely n — 2 photon’s transverse 
masses set at the fundamental mode, Mj = fee, * = 2 ~ n — 1. The result is expressed by the inequality, 

^9 

We may define the largest McQ3 trigger energy Wc by taking n = 3,uj = kc, 


Uc = 


£eg 

2 




(13) 


The trigger effective mass Mq is often taken to be fee in the present work. 
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On the other hand, a massive neutrino of mass m satisfies the dispersion relation of the energy and the 
momentum, E = \Jp^ + A big difference is in transverse components: the neutrino transverse 

momentum may vanish: pj_ = 0. This gives an allowed region for RENP; 


Wo < 


ceg 


— Mq) . 


RENP thresholds of neutrino pair production of mass m* start at 


/ _ 
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(14) 


(15) 


This location is shifted to the higher energy side by an amount MQ/(2eeg) from the free space value ua. 
Mq = kc for the trigger fundamental mode. 

In both of inequalities, eq.([T2|) and ea. IfTTIl . the special cases of equalities give the boundary curves of 
two allowed regions. Thus, there exists a gap between RENP allowed regions and McQn allowed regions, if 


(n — l)kc > 2mo , n > 3, 


(16) 


where mo is the smallest neutrino mass. The extra kc factor here arises from the smallest energy of u equal 
to kc- In the gap region RENP is QED background-free. The argument breaks down at n = 2, where the 
usual relation w = wq = eeg/2 + k^ji^ecg') holds [T^ . 

Thus, all QED backgrounds are rejected once it is rejected at the smallest n value equal to 3. We 
illustrated the background-free region in FigQ. 


RENP and McQ3,McQ5 constraints 



Figure 4: QED background-free regions of Xe RENP in (wo,w) plane. Two curves give upper boundaries 
of McQn (n=3 in solid black, =5 in dashed red) allowed regions. Boundaries are given for the fundamental 
mode (nf + = 1) and those for excited modes (not shown here) are all within these allowed regions. 

The background-free RENP is possible in the indicated horizontal wq— region in red up to the arrow point. 
Assumed parameters are teg =8.315 eV, and the size d = lO^um of photonic-crystal fiber. 


RENP rates 
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We consider RENP from 6s^Pi (8.437 eV) [Hj. We need data of 7s^Pi (10.593 eV), and its A-coefRcient, 
A{7s — 7 > 5p) = 8.51 X 10^s“^, since the Coulomb excitation bridging between the nucleus and the valence 
electron requires another state of the same quantum number of Gs^Pi. See |10] for more details. Xe 
RENP rate from the nuclear pair emission, both in free space and in host guides (assuming the trigger field 
to be an allowed propagating mode), is given by 


t) = TQFx{uj)I{u:)r]^{t ), 


Pq ~ 54mHz—^( 
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V , lOOMHz 


eV^lO^^cm"^ lO^cm^ eV^ 
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(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


^7s ^7>p \/37r {(-7s (6s 4“ ^)(f^6s (5p 

where 6 m = 1 for the Majorana neutrino and 6m = 0 for the Dirac neutrino [1]. 

Both in wave guides and photonic crystals, overlapping integrals containing mode functions need to be 
evaluated. They typically take the form. 


1= [ dxcos^^e-*^"^ 
Jn d 


riyTry 


dy sin “ e 
a 


( 22 ) 


for TE modes. Here Pj’s are transverse components of momentum of the neutrino pair. Let us restrict the 
trigger to the fundamental TEqi mode: ni = 0,77-2 = 1. An elementary calculation gives 


|j|2 ^ ^g^ sin(P,,d/2) ^ 2 ^Py cos{Pyd/2) 


P. 


Px,y — Px.v 4 " P: 


x,y ' ^x,y • 


(7r/d)2 - P2 ^ 

In the large size limit of d —>• oo, the functions here are all reduced to the delta functions, 


(23) 


|/|2 ^ 2{277)‘^dH{Px) [6{Py - ^) + 6{Py + ^)) . (24) 

The back to back momentum balance of transverse momenta {p±)i of the neutrino pair is broken by an 
amount, 4=7r?7o/d according to this formula. Thus, the spectrum rate calculated from these delta functions 
coincides with the continuum limit calculation done as if the trigger were slightly off-axis. Since the squared 
matrix element of RENP is a smooth of the transverse momentum around P_j_ = 0, the approximation of 
large size limit, d —oo, is justified. 

As is evident in the contour map for the background-free region, all cuo—region of uq > ojc (<^c defined by 
m ) for RENP belongs to the McQn (n > 3) background-free region. The background-free RENP spectrum 
consists of a portion of the RENP spectrum truncated at ujc < uiq < ujqq, with ujqq determined by ea. ffT^ 
taking the smallest neutrino mass mo- The full spectrum in the background-free region is shown in Fig([5]). 

Finally we show in Fig([6]) RENP spectrum in McQn (n > 2) background-free region, assuming lO/im 
size of photonic-crystal fiber, to show all three thresholds uJii,i = 1,2,3 are measurable. Note that the 
experimental resolution of wq is made possible by the excellent frequency resolution of lasers, typically of 
order 10“^ eV or 100 MHz. 













Xe RENP spectrum:Majorana NH vs IH 
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Figure 5: Xe spectral rate in the background-free region for the Majorana cases of the smallest mass, mo = 
1 and 5 meV, and NH (in solid black and in dash-dotted blue) and IH (in dashed red and in dotted black) 
. Xe gas density of 1 x 10^*^cm“^, the target volume Icm^, and rj^i = 10“^ are assumed. The assumed size 
of a host guide is 100 fim. 


Xe RENP spectrum:Majorana NH vs IH 
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Figure 6: Xe spectral rate in the background-free region for the Majorana cases of smallest masses, 1 and 5 
meV, and NH (in solid black and in dash-dotted blue) and IH (in dashed red and in dotted black). Xe gas 
density of 1 x 10^°cm~^, the target volume Icm^, and = 10“^. The assumed size of a host guide is taken 
as 10 /im. Three neutrino pair emission thresholds are clearly visible. 
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In summary, we demonstrated that QED background-free RENP may be realized in wave guide or in 
some type of photonic crystals. With a large enough cutoff of transverse wave number three neutrino pair 
emission thresholds are detectable, using the largest nuclear monopole contribution. 

We should like to thank J. Arafune for raising the QED background problem in the correct way. 


4 Appendix: Propagating field modes in wave guides and photonic crys¬ 
tals 


Allowed modes in rectangular wave guides are given by discretized transverse stationary waves. Eor example, 
the simplest TE mode propagating in z-direction have electric field components of the form [6], 


^ UajTTX n^vru , 1 . „ 

Ex oc cos-sm 


ocsin^^cos^e'^'^A 

y „ b 


(25) 

a b " a b 

with Ez = 0. Ux^riy are both integers. With a < b the lowest fundamental modes are TEqi and TMqi of 
the eigenvalue — A 2 = —d’^jdx^ — /dy^ is {t: jEf'. Eor simplicity of notations we shall take the limit of a 
square wave guide, a = b = d, however without further consideration of degeneracy of modes. These mode 
functions are orthogonal, and form a complete set for the general expansion of propagating light waves. The 
lowest excitation mode of TEio is given by a non-vanishing held of 

Ejj . ,7r , ji,_ ^ „ //TT, 

■ sm(—7/je 


Ey{x,y,z-,t) = 


^ikzZ 


Ex{x,y,z]t) = 0, 




(26) 


with V = d^L the volume of the target. The normalization of helds is determined by the condition that the 
energy density of this held per unit length is u)/L. 

Allowed projected momentum k_\_{= '^{nx-,ny)/d for TE modes) of the light held forms on a lattice given 
by discrete points in the projected two dimensional space orthogonal to the trigger axis, denoted in general 
by nn/d. Unless the projected momentum falls on these lattice points, the photon emission is forbidden 
in wave guide. We may introduce an effective cutoff wave number by kc = vr/d: transverse momenta of 
|/cj_| < kc are forbidden in the wave guide. The important point for the background rejection is that this 
cutoff acts as an effective mass to parallel momentum component k\\ of w = \Jk‘^\+k‘^^. 


Eigenvalue problem in photonic crystals 

We shall recapitulate basic features of eigen-modes that propagate in cylindrical hbers of photonic 
crystals. 

Eirst, the fundamental eigenvalue equation is set up for the magnetic held [7j, 

oj[k)e{r) 

{ik + V)x -h V) X ttg(f) = a;2(fe)ttg(r) , (28) 

where e(r) is the dielectric function which may depend on spatial places. The linear operator acting the 
vector function is hermitian. 
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In the case of cylindrical photonic-crystal fibers, one may use the cylindrical symmetry of the dielectric 
function e(r), to factor out the cyclic coordinate dependence, 

= p=^x^ + y\ (29) 

where kz is an arbitrary real number and n is an integer. We shall not write down an ordinary differential 
equation for the vector function hk^^n{p)- This is the eigenvalue problem in one dimension (ID) due to the 
cylindrical symmetry. 

When the eigenvalue problem is applied to a multi-layer dielectrics of alternating dielectric constants 
with a large-size air hole in the center [8], there exist band gaps and their dispersion relation u}{kz) provides 
the effective cutoff of order n/d, where d is the cladding periodicity, the hole radius being taken comparable 
to d. Some details are shown in [9]. 

5 Appendix: Spin current contribution from ^P 2 (8.315 eV) de-excitation 

For completeness we show in Fig([7]) the spin current contribution using the formulas of [11] . Although rates 
are much smaller than the nuclear monopole contribution of preceding figures, this case has a potentiality 
to access the Majorana CPV phase determination. 
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